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^ . ABSTRACT 

| Observational manifestations of possible breaking of spin-statistics relation for 

£S| ■ neutrinos are considered. It is argued that bosonic neutrinos may form cosmolog- 

; I | ical cold dark matter, improve agreement of BBN predictions with observations, 

Q-i' make operative Z-burst model of ultra-high energy cosmic rays, etc. Restrictions 

| for an admixture of bosonic component to neutrino which follow from double 

\q ■ beta decay are discussed. 
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1. Introduction 

> . 

There is an impressive symbiosis of the "weakest" and lightest of the known mas- 
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sive particles, neutrino, and Cosmos. Cosmology and astrophysics allow to study 



■"sj- \ neutrino properties with an accuracy which, in many cases, is unaccessible in direct 



terrestrial experiments an d, vi ce versa, neutrino helps to resolve some cosmological 
and astrophysical mysteries'!***'. 

Cosmology allows to put a very stringent upper bound on neutrino mass at the 
level of about 1 eV, for a review see e.g. papers"*. The bound on the amplitude of 
Ph. possible right handed currents and the mass of right intermediate bosons found from 

the analysis of Big Bang Nucleosynthesis (BBN) is orders of magnitude better than 
those obtained in laboratories. The same is true for neutrino magnetic moments and 
mixing of the usual active neutrinos with hypothetical sterile ones which are restricted 
by BB]>P and stellar evolution®. 

On the other hand, neutrino contributes to cosmology providing hot dark matter, 
but not the necessary cold one, if physics is normal. Neutrino may be related to dark 
energyE3, anc [ b e responsible, at least partly, for ultra-high energy cosmic rays beyond 
GZK cut-off through the Z-burst mechanisnJ-3 . The large mixing angle solution to 
the solar neutrino anomaly excludes noticeable lepton asymmetry of the universe^*. 

These lists are far from being complete but this is not the main subject of this 
talk. Instead of these rather well known topics I would like to talk about new, though 
quite speculative issues, related to effects of possible breakingof neutrino statistics 
in cosmology. The content of the talk is strongly based on ref™. 

The first question is why neutrino? First of all, neutrino is the only known particle 
indicating to new physics. As is known from observation of neutrino flavor oscillations 
leptonic flavor charges, electronic, muonic, and tauonic, are not conserved. Neutrino 
is the only observed particle for which Majorana mass is possible and, as a result 



total leptonic charge could be non-conserved. There are experimental indications to 
a possible leptonic charge non-conservation from neutrinoless double beta decaj^. 
It may mean that neutrino communicates with a hidden sector of the particle world 
and is a messenger of new physics from the hidden sector. It could be that no sacred 
principles are respected in the hidden world and neutrinos brings us exotic possibilities 
of breaking CPT theorerrP^ or Lorentz invariance, which are actively discussed in 
the recent years. Since at the present days cosmology is quickly becoming precise 
science, maybe cosmos will bring through neutrinos new surprising physics. 

The most exciting possibility which, in particular, may lead to violation of CPT 
and Lorentz invariance and even to much more drastic consequences is a breaking of 
the spin-statistics relation for neutrinos. Ironically the particle brought to this world 
by Pauli may violate Pauli exclusion principle. In fact Pauli and Fermi repeatedly 
asked the question if spin-statistics relation could be not exact and electrons were a 
little bit different. 

Possible violation of exclusion principle for the usual matter, i.e. for electrons and 
nucleons was discussed in a number of papers at the end of the 80tfP^. Efforts to find 
a more general than pure Fermi-Dirac or Bose-Einstein statistics^ 3 were taken but 
no satisfactory theoretical frameworks had been found. Experimental searches of the 
Pauli principle violation for electrons and nucleons have also given negative 
results. 

If one assumes that spin-statistics relation is broken while otherwise remaining 
in the frameworks of the traditional quantum field theory then immediately several 
deep theoretical problems would emerge: 

1) non-locality; 

2) faster-than-light signals; 

3) non-positive energy density and possibly unstable vacuum; 

4) maybe breaking of unitarity; 

5) broken CPT and Lorentz invariance (as mentioned above). 

Either these consequences (if they indeed were realized) exclude any violation of 
spin-statistics theorem and discussion of this violation should be forbidden or they 
open an exciting space for further research and development. An answer to that is 
first of all a matter of experiment which may either exclude or confirm the drastic 
assumption of breaking the spin-statistics relation. As for observational manifesta- 
tions of the mentioned phenomena they should be weak because they are induced by 
weakly interacting neutrinos and, moreover, in higher orders of perturbation theory. 

Perturbative expansion of the scattering matrix has the well known form: 



where T{...} means time-ordered product of operators inside brackets. Lorentz in- 
variance is ensured if Hamiltonian, Tl, is bosonic operator, i.e. it commutes with itself 




(1) 



if separated by the light cone, see e.gr^. However, for bosonic or partly bosonic neu- 
trinos the effective Hamiltonian responsible even for the simple reaction e+p <-> n + u 
is not bosonic and observables do not commute for space-like separation and locality 
breaks. Presumably unitarity is maintained because Hamiltonian remains Hermitian. 

Another possibility that Hamiltonian/Lagrangian approach and least action prin- 
ciple are applicable only approximately and theory is drastically modified, while the 
observable effects may still be small. 

So let us postpone discussion of (non-existing) theory and consider phenomenology 
of neutrinos obeying Bose or mixed statistic J^. What can we buy for this price? 



2. Dark matter . 



It is well known that the usual fermionic neutrinos cannot form cosmological cold 
dark matter for any spectrum of primordial fluctuations and arbitrary self-interaction. 
This conclusion is based on the Tremain-Gunn bouncP^ which does not allow to 
fill galaxies with sufficiently many light fermions (satisfying the Gerstein-Zeldovich 
boun to account for the observed hidden mass. Thus we face the following dark 
matter dilemma: 

1) new particles and old (normal) physics 

2) old particles (neutrinos) and very new physics. 

To make the cosmological cold dark matter neutrinos must form Bose condensate 
in the early universe. To this end a very large lepton asymmetry is necessary with 



n„ - n,-. 



100 (2) 



Such asymmetry might be created in a version of the Affleck and Din<Pl scenario. 

A large asymmetry allows to fill up the present day universe by a huge number of 
cosmic neutrinos such that they would be able to make all CDM, Qcdm ~ 0.25, if 

n u ~ 10 4 cm" 3 (3) 

The spectrum of cosmic background neutrinos, if they are bosonic, would be very 
much different from the fermionic ones because the Bose equilibrium distribution has 
the form: 

U b = r , F \ /T : +C5(k), (4) 

exp[(E - ii v )IT - 1 

where /i = m v is the maximum value of chemical potential of bosonic neutrinos equal 
to their mass. The condensate amplitude C does not depend upon neutrino energy 
but may depend upon time. One sees that the bulk of the condensed neutrinos is 
cold. In galaxies the neutrino number density would be about 6 orders of magnitude 
larger than the average cosmological number density, i.e. 

n {gai) „ 10 io ( mjy /o.i eV) cm" 3 . (5) 



Structure formation with Bose condensed light bosons with the usual integer spin 
was considered in refff^. The model well reproduces the essential features of the 
observed large scale structure. Since the picture is spin independent the same must 
be true as well for bosonic neutrinos. 

The results and numerical estimates presented in this section are true for purely 
bosonic neutrinos however, as we see below, experiments on two neutrino double beta 
decay seem to exclude 100% breaking of statistics and at least some fermionic fraction 
must be present in a neutrino. It makes the model noticeably more cumbersome, but 
less vulnerable. 

3. Equilibrium distribution for mixed statistics 

The statistics dependent term in kinetic equation for the reaction 1 + 2 <-> 3 + 4 
has the form 



F = /i(pi)/ 2 (p 2 )[l ± /s(p 3 )][l ± U{pd) ~ fs{ps)Mp^ ± /i(Pi)][1 ± / 2 (p 2 )] (6) 



where fj is the distribution function of particle j. This expression is valid in the 
case of T-invariant theory when the detailed balance condition is fulfilled. Since T- 
invariance is broken, kinetic equation is modified but the equilibrium distributions 
fj eq ^ remain the same canonical Bose and Fermi ones as in T-invariant theory 
This statement is based on the unitarity of S'-matrix. If the spin-statistics relation 
is broken, as a result the unitarity may also be broken. If this is the case, then a 
breaking of T-invariance may create large deviations from the standard equilibrium 
distribution functions of neutrinos. 

In what follows we neglect complications related to a violation of T-invariance. In 
the case that neutrino obeys pure Bose statistics its equilibrium distribution is given 
by the standard Bose form (J3J). Indeed, it is easy to see that for this distribution F 
vanishes and together with it the collision integral vanishes too. However, the form 
of equilibrium distribution for mixed statistics is not so evident. It depends upon an 
assumption about F for particles obeying mixed statistics. We do not have rigorous 
arguments in favor of one or other form for F and as a reasonable guess assume that 
the factor depending upon the neutrino statistics in F changes as 



where c = cos 7, s = sin 7 and 7 is some mixing angle characterizing admixture of 
wrong statistics. 

Another possibility for description of mixed statistics in kinetic equation could be 



However, it is easy to see that these two seemingly reasonable possibilities (jjj) and 
(jHJ are identically equivalent. In both cases 




(7) 






(9) 



where 

K = c 2 - s 2 = cos 27 (10) 

We call k the Fermi-Bose mixing paramete One can check that in the case 

of mixed statistics introduced to kinetic equation according to Q the equilibrium 
distribution takes the forrr!^: 

/<«*> = [exp(E/T) + K }- 1 . (11) 

where k runs from +1 to —1 corresponding respectively to Fermi and Bose limits. 
The intermediate value k = corresponds to Boltzmann statistics. 

If — 1 < k < 0, the maximum value of the chemical potential may be bigger than 
the neutrino mass: 

/1 (mai) = m,-Tln(-K) (12) 

Bose condensation might take place for negative k only. 

Another possible form of a modification of statistics dependent factor in kinetic 
equation would emerge if we assume that there are two neutrino fields with the same 
mass and different statistics, fermionic and bosonic. The Lagrangian would always 
depend upon two independent field operators in the combination: 

ipv = cip b + sijjf, (13) 

where ipj^ are respectively bosonic and fermionic operators. In this case kinetic 
equation would contain two different terms: 

c 2 ff(l-ff) (14) 

and 

s 2 m - h) (15) 

Equilibrium distributions would be canonical ones, e.g. for vanishing chemical poten- 
tial they are: 

f fJb = 1/ [exp (E/T) ± 1] (16) 

but the number of states in equilibrium becomes doubled. On the other hand, the 
probability of e.g. neutron beta-decay remains the same as in the standard theory 
because c 2 + s 2 = 1. 

4. Big Bang Nucleosynthesis 

The impact of purely bosonic neutrinos on BBN was considered earlier in paper^l. 
The effects of mixed statistics described by the equilibrium distribution (|TT|) were 
calculated in our wor The equilibrium energy density of bosonic neutrinos at 
T ^> m u is 8/7 of the energy density of fermionic neutrinos and thus the change of 



statistics would lead to an increase of the effective number of neutrino species at BBN 
by AN U = 3/7 (for three neutrinos). On the other hand, a larger magnitude of the 
neutrino distribution function and the fact that it enters kinetic equation (see (JBJ)) as 
(l + f u ) instead of (1 — f v ) makes the weak reactions of neutron-proton transformations 
faster and hence the n/p freezing temperature becomes lower. This effect dominates 
and as a result the effective number of massless species becomes smaller than 3. 
According to the calculations of refP^ the effective number of neutrino species in the 
case of pure Bose statistics becomes N e ff = 2.43, practically independently on the 
value of the baryon-to-photon ratio rj = ris/n^. 

The effective number of neutrino species determined by the comparison of the 
calculated abundance of primordial 4 if e with the standard result is presented in the 
upper panel of fig. [T] as a function of k. However, the effect of change of statistics 
cannot be described by a simple change in N v if other light elements are included. In 
the lower panel of fig. 1 the relative changes of the abundances of 2 H, 4 He, and 7 Li 
are presented. As expected the mass fraction of A He drops down, while the amount 
of 2 if goes up. A higher deuterium abundance can be explained by a slower conver- 
sion of deuterium to heavier elements due to fewer neutrons and faster cosmological 
expansion at T rs 0.8 * 10 9 when the light elements have been formed. 

At k = — 1 we find for 4 He: Y p = 0.240, which makes much better agreement with 
the value extracted from observations (for a review of the latter see e.g.^S). Different 
helium observations yield different results, e.g., ref. finds Y = 0.238 ± 0.002, and 
ref. finds Y = 0.2421 ± 0.0021 (la, only statistical error-bars). These results are 
shown in figure 2 as the skew hatched (yellow) region. Whether the existing helium 
observations are accurate or slightly systematically shifted will be tested with future 
CMB observations^. 

The amount of 2 H rises at most to Xi H fXn = 2.5 • 10~ 5 and the agreement 
between BBN and WMAP data remains good, bearing in mind the observational 
uncertainties. Primordial 7 Li drops down to Xi^jXn = 4.55 ■ 10~ 10 , again slightly 
diminishing the disagreement between theory and observations. 

We see that at the present time BBN does not exclude even a pure bosonic nature 
of all three neutrinos. Furthermore, the agreement between the value of the baryonic 
mass density, rj, inferred from CMBR and the predicted abundances of 4 He, 2 H, and 
7 Li becomes even better. In other words, in the standard BBN model there is an 
indication of disagreement between observations of 4 He and 2 H - they correspond to 
different values of rj with the observed abundances of 4 He indicating a smaller value 
than given by CMBR, while 2 H agrees with CMBR. Motivated by these results the 
value of AN U = —0.7 ± 0.35 was suggested in ref. ^lL ) . In the case of predominantly 
bosonic neutrinos, as discussed above, the discrepancy between 2 H, 4 He, and CMBR 
disappears. 

The results presented in this section are obtained for negligible chemical poten- 
tial of electronic neutrinos. On the other hand, formation of cosmological neutrino 




Figure 1: Upper panel: the change in the effective number of degrees of freedom which corresponds 
to the change of the 4 He abundance as a function of the effective Fermi-Bose parameter k. Lower 
panel: the relative change of the primordial abundances of deuterium, helium-4, and lithium-7, as 
functions of k. We take 77 = ns/n 7 = 6.5 ■ 1CP 10 . 




Figure 2: Upper panel: the ratios of abundances of different elements in the cases of purely bosonic 
neutrinos with respect to the standard fermionic case as functions of the baryon number density, 
rj. The vertically hatched (cyan) region shows the WMAP 2a determination of rj. Lower panel: 
the absolute abundance of 4 He as a function of rj for the purely bosonic, Boltzmann, and fermionic 
neutrino distributions, corresponding to k = — 1,0, +1 respectively. The two skew hatched regions 
show the observation of primordial helium from ref. 24 ) (lower, yellow) and ref. 25 -* (upper, magenta), 
which marginally overlap at la. 



condensate discussed in sec. 2 demands the maximum value of \i v given by eq. (fT2"|) . 
As is knowrP, BBN allows at most fi/T = 0.07 for any neutrino flavor. This implies 
n > 0.9. Such a large admixture of a wrong bosonic state to v e is most probably ex- 
cluded by the data on double beta decay (see below, sec. 7). However, one may still 
hope to save the neutrino cold dark matter (z/CDM) if the mixing angle determined 
from the decay is different from that that enters neutrino kinetics (see discussion in 
sec. 8. Another possibility to save z/CDM is to assume that the chemical potentials 
of Up or \x T are much larger than that of v e . The latter are very weakly restricted by 
BBN and only the large mixing between neutrino flavors equalizes all chemical po- 
tentials. However, the change of neutrino statistics may lead to a different refraction 
index in the primeval plasma and to suppression of the transformation of z/ M T into u e . 
There is also a more conventional way to suppress neutrino flavor oscillations in the 
primeval plasma introducing neutrino coupling to light pseudo-goldstone boson, Ma- 
joron. The effective potential of neutrinos induced by the Majoron exchange would 
suppress flavor transformations in the cosmological plasmsP^J. This would allow to 
have large chemical potentials of v^ T and small chemical potential of v e . 

5. Astrophysical consequences 

Neutrino statistics plays key role in the environments where neutrinos form dense 
degenerate gases. Direct test of the "bosonic" nature of neutrinos can be provided by 
precise measurements of the neutrino energy spectrum from supernova. Generically, 
the spectrum of bosonic neutrinos should be more narrow. To establish the difference 
one needs to measure the spectrum both in the low, E < 3T, and in the high, E > 3T 
energy parts. 

A violation of the exclusion principle can influence dynamics of the SN collapse. 
According to the usual scenario at the initial stages (formation of the hot proto- 
neutron star) the neutronization leads to production of high concentration of the 
electron neutrinos which are trapped in the core. The chemical potential of these 
neutrinos (due to the Pauli principle) can reach 70 - 100 MeV. These neutrinos heat 
the medium and diffuse from the core. Violation of the Pauli principle allows for the 
neutronization neutrinos to be produced with lower energies. These neutrinos escape 
easier the star leading to faster cooling and lower central temperatures. The evolution 
of the lepton number would change as well. 

High neutrino density in the condensate (especially if an additional clustering 
occurs) enhances the rate of the Z°-bursts produced by the annihilation of the ultra 
high energy (UHE) cosmic neutrinos on the relic neutrinos This in turn, enhances 
production of the UHE cosmic rays, and may help to explain the cosmic ray evens 
above the GZK cut-off. 

Charge asymmetric neutrino condensate may produce a strong refraction of the 
high energy neutrinos from remote sources (active galactic nuclei, gamma ray bursters). 



Apart from lensing, one may expect a substantial impact on neutrino oscillations^*. 

Since the density of dark matter in galaxies is about 6 orders of magnitude larger 
than their average cosmological energy density, a condensation of cold neutrinos 
around the Earth might have an effect on the end point of the beta decay spectra, in 
particular, in the tritium decay experiments on search for neutrino mass . 



6. Double beta decay 

In contrast to electrons and nucleons which form atoms and nuclei, where the 
effects of statistics are of primary importance, it is difficult to observe processes with 
identical neutrinos. A realistic reaction for the test of neutrino statistics can be the 
two-neutrino double beta decay, 

A -> A' + 2v + 2e~ (17) 

(or similar with production of antineutrinos and positrons). The probability of the 
decay as well as the energy spectrum and angular distribution of electrons should be 
affected by the change of neutrino statistics. 

To have a formalism for description of identical neutrinos one needs to specify 
operators of neutrino creation/annihilation. We assume that they consist of two parts, 
fermionic, /, and bosonic b for operators of annihilation, a — f + b. Its Hermitian 
conjugate could naturally be the operator of neutrino creation. Correspondingly we 
define one neutrino state as: 

\u) = a+|0) = Ci ;+|0> + Sl S+|0) = c\f) + s\b) (18) 

where |/) and \b) are respectively one particle fermionic and bosonic states and c\ = 
cos S and si = sin 5. It would be natural to expect that 5 is equal to 7 introduced 
above in eq. (JJJ) but we cannot prove it formally. 

To describe the two-neutrino state one needs to specify the relevant commutators 
which, we postulate, have the following form: 

fb = e%f, f + b + = e% + f + , fb + = e-% + f, f+b = e~%f + , (19) 

where <fi is an arbitrary phase. The two neutrino state is natural to define as 

\h,k 2 ) = a+a+\0} (20) 

The matrix element of the decay of nucleus A into 2u + 2e + A' may be possibly 
taken in the usual way: 



A 2 p = (h, h,2e, A 1 J (Tx X (rx2^v (^1)^2(^2) M(xi,x 2 ) 0, A). (21) 
After making the necessary commutating according to eq. (|T9*|) we obtain: 



A 2P = A. 



4 + c(s( (1 -cos0) +A + s? + cfsf(l + cos0) . (22) 



where v4_ and A + are respectively antisymmetric (fermionic) and symmetric (bosonic) 
parts of two neutrino emission. It is easy to see that the amplitude can be parametrized 
as 

A w = cos 2 x A. + sin 2 X A+, (23) 

where cos 2 x = c 4 + c 2 s 2 (1 — cos <fi) and sin 2 x — s \ + c \ s \ (1 + cos <f>). The probability 
of the double beta decay integrated over neutrino momenta evidently does not contain 
interference between A + and A_ and is equal to: 

W tot = cos 4 X W. + sin 4 X W + , (24) 

where W± are proportional to |A±| 2 . 

The probability of decay into unusual bosonic neutrinos is proportional to the 
bi-linear combinations of the type K m K n , K m L n , L m L n , where 

Km = [Em — Ei + E el + E u iY x — [E m — Ei + E e2 + E v2 \~ x , 
L b m = [Em — Ei + E e2 + E u i] 1 — [E rn — Ei + E e \ + E v ^~ l . (25) 

Here the upper index b indicates that the results are applicable to bosonic neutrinos, 
Ei is the energy of the initial nuclei, E m is the energy of the intermediate nucleus 
state m, E ei , and E ui are the energies of electrons and neutrinos respectively. The 
minus signs between the two terms in the above expressions are due to the bosonic 
character of neutrinos; in the case of fermionic neutrinos we would have plus signs^^. 
For electrons we assume the normal Fermi statistics. 

In the case of + — > + transitions the combinations K m and L m can be approxi- 
mated by 

K b E e2 — E el + E u2 — E vl b E el — E e2 + E v2 — E vX . , 

m ^ (E m -E i + E /2y ' m ^ (E m -E l + E /2Y ' 1 j 

whereas for the fermionic neutrinos 

9 

»1 m 771 771 , 771 In' \ / 

Here E /2 = E e + E v is the average energy of the leptonic pair. Appearance of the 
differences of the electron and neutrino energies in eq. leads to a suppression 

of the total probability. It also modifies the energy distributions of electrons. The 
probabilities of the transitions + — > 2 + are proportional to the combinations (K m — 
L m )(K n - L n ), where 

(XC* Tb \ ^ 2(Ee2 - E e i) 

{Km ~ LJ ~ (E m -^ + Eo/2) 2 - (28) 

In the case of fermionic neutrinos the combination has an additional factor (E u2 — 
E u i)/(E m — Ei + Eq/2) and the suppression is stronger. 



A simple estimate shows that the probability of + — > + -transition for bosonic 
neutrinos is suppressed by 1/250 for 56 Ge and by 1/10 for 100 Mo. Theoretically the 
total decay rate is known with the accuracy within a factor of few. This probably 
allows to exclude a 100% bosonic neutrino. However, the fraction of bosonic neutrino 
can still be very high. According to our preliminary calculation^ 2 ' 2 ^' the value of the 
mixing angle can be as large as: 

sin 2 X < 0.7 (29) 

For + — > 2 + the situation is opposite: bosonic neutrinos are more efficiently 
produced. However, no interesting bound is obtained in this case because the statistics 
for these decays is much lower. 

One can use the data on the spectrum of the emitted electrons, either single 
electron spectrum or distribution over the total energy of both electrons. The spectra 
do not have any noticeable ambiguity related to unknown nuclear matrix elements 
and the present day accuracy is at the level of 10%. Potentially their analysis may 
improve the above quoted limit (|29jl or indicate the existence of a "bosonic" admixture 
to neutrinos. Some already observed anomalies may be interpreted as hints supporting 
the latter. 

Unfortunately we cannot say at the present stage how the Fermi-Bose parameter 
introduced above (fTUJ) is related to the mixing angle \. Even if we assume that the 
mixing angle in neutrino kinetics (J7J) is the same as in the definition of neutrino states 
(11811 . the unknown value of the angle <fi which enters the commutation relations (|19p 
and upon which depends the angle x (pH| makes the relation between k and \ rather 
arbitrary. 



7. Theoretical problems and discussion 



Mentioned above ambiguities are related to intrinsic problems of formulation a 
theory with mixed statistics. Working at a naive level, as we did above, it is even 
difficult to define the properly normalized particle number operator. According to 
eqs. (jTSJl and (}2Tjj) it is natural to define the n identical neutrino state as 



n) 



i) = {cj + + s 1 b + y io> 

The normalization of this state is 



(n\n) 



2(n-l) 

6 i 



n\s\ + (n- l)\c\ 



( sin (n<j)/2) 



\ sin (0/2) / 

If we introduce the particle number operator in the usual way: 



(30) 



(31) 



h = a a, 



(32) 



then its diagonal matrix elements would be 



2 o <^n — 1) sinn0/2 



nn\s l + 2n\c 1 s 1 cos : — — - + 

2 sin 0/2 



/ I " I \ 2(n-l) 

4H + (-i)'«-«)(^gi) 



(33) 



The particle number operator, as introduced above, has reasonable and self-consistent 
interpretation only for the case of pure statistics, while for mixed statistics it even 
does not commute with Hamiltonian if the latter, or operator n, or both are not 
somehow modified. 

There are no problems with reactions where only one neutrino is involved, but 
serious difficulties may arise with two neutrino reactions, as e.g. ul — > vl or vv — > //, 
even if the participating neutrinos are not in an identical state. The amplitude of ve- 
elastic scattering in the usual approach is given by the expansion of the T-exponent 
of the action and is described by two diagrams differing by an interchange of emission 
and absorption points. If taken literally, the diagrams with W /± -exchange would give 
vanishing amplitude for purely bosonic neutrinos a . In this case only Z-exchange 
would contribute to i/ e e-scattering and the cross-sections of v^e and z/ e e-scattering 
would be equal. Reactor neutrino experiments are consistent with the standard value 
of v e e cross-section and seem to exclude the possibility of purely bosonic v e . Using 
these data one can put a rather strong bound on bosonic admixture to electronic 
neutrino. On the other hand, perturbation theory with non-bosonic Hamiltonian 
may need to be modified and the above conclusion of vanishing of the amplitude of 
scattering of pure bosonic neutrinos on electrons would be invalid. 

It is unclear if all these problems can be resolved in a simple way or drastic 
modifications of the underlying theory is necessary, which is a nontrivial task because 
the observed consequences of the theory must not be destroyed. 

The presentation in the previous sections and in the original papeJ^ was on pure 
(and poor) phenomenological level. For example if neutrinos have mixed statistics 
then in double beta decay the symmetry of the final state of neutrinos is mixed: 
symmetric with wight a + and antisymmetric with the weight a_. It seems plausible 
that these weights are respectively cos 4 \ an d sin 4 \ as argued in the previous section, 
eq. (|24)1. simply on the basis on the normalization arguments. 

Similar reasoning is possible for kinetic effects, eqs. ()7I8|) . There are no rigorous 
theoretical arguments in favor of such description but the result for the equilib- 
rium distribution in the case of mixed statistics looks quite beautiful. Moreover, the 
fact that two "reasonable" (or natural) ways of description (J7J) and (jHJ) give the same 



a This was noticed also by F. Vissani at this Conference. 



result is an argument in favor of their validity. 
8. Conclusion 

There is no consistent theoretical frameworks for description of mixed neutrino 
statistics and even in the case of purely bosonic neutrinos the fermionic property of 
the Hamiltonian would make a possible future theory quite unusual if it will ever be 
formulated. Still independently on theory there could be some predictions testable 
by experiment. So to summarize we will conclude that: 

1. The suggestion of bosonic or mixed statistics for particles (neutrinos) with half 
integer spin looks exciting but opens a Pandora box of serious theoretical prob- 
lems, which may be impossible to resolve without revolutionary modification of 
the standard theory. Such modification looks especially difficult in the case of 
mixed statistics. 

2. The suggested mixture of statistics allows to break plenty of sacred principles, as 
e.g. Lorentz invariance, CPT-theorem, locality, etc, which are actively discussed 
now. 

3. Bosonic neutrinos open a possibility of making all cosmo logical dark matter out 
of neutrinos in accordance with Occam's razor: "Plurality should not be posited 
without necessity." 

4. "Bosonization" of neutrinos leads to effective number of neutrino species at 
BBN smaller than 3 and makes an agreement of the BBN calculations with the 
data noticeably better. 

5. Analysis of accumulated and accumulating data on two neutrino double beta 
decay could restrict the admixture of wrong statistics to neutrinos or to indicate 
a violation of spin-statistics relation. 

6. Last, but not the least, if the validity of spin-statistics theorem has been stud- 
ied for the usual matter, electrons and nucleons it surely worth studying for 
neutrinos. The possibility that statistics is modified for neutrinos seems more 
plausible because neutrino is a natural particle to be a messenger from hidden 
sector of physics where some principles respected in our world can be violated. 
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